Gaussian Schell Source as Model for Slit-Collimated Atomic and Molecular Beams 
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We derive the spatial coherence and intensity profiles of beams emerging from two consecutive 
collimating apertures, and compare our results with data. We show how to make a Gaussian Schell- 
model (GSM) beam by assuming Gaussian apertures. Then we compare the intensity profile, the 
transverse coherence width and the divergence angle of a GSM beam with those same properties 
of a beam that is collimated with two hard-edged slits. The GSM beam formulae are simpler, and 
offer an intuitive way to understand how partially coherent beams interact with interferometers. 



I. INTRODUCTION 

Perhaps the simplest way to form a beam of particles 
is to place two consecutive apertures in front of a source. 
As illustrated in Figured! a first aperture placed directly 
in front of an extended incoherent source helps define the 
source. Then a second aperture, farther away, collimates 
the beam. In many experiments with X-rays, gamma 
rays, atoms, molecules, neutrons or electrons a simple 
ray picture can adequately predict the intensity profile 
of such beams. But more care is needed to model the 
spatial coherence that is so crucial for understanding the 
behavior of interferometers. 

A Gaussian Schell-model (GSM) beam is a particular 
type of partially coherent wave field that is often used to 
model beam coherence properties (transverse coherence 
lengths and wavefront radius of curvature) as well as in- 
tensity profiles. In fact, we have learned quite a bit about 
three-grating interferometers with GSM beams [1]. The 
major benefit of GSM beams is that the Huygens-Fresnel 
propagation can be described analytically in many cases. 
However, the question has been raised, Ts the GSM a 
realistic way to model the beams in matter wave inter- 
ferometers?' The purpose of this paper is to answer this 
question in the affirmative. We show how to make a 
GSM beam and then we describe how the GSM beam's 
properties compare to a slit-collimated beam made with 
hard-edged slits. 

Furthermore, we show how a simple Gaussian Schell- 
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Figure 2: Actual intensity profile of a beam emanating from 
collimating hard-edged slits (circles) and Gaussian apertures 
(squares), compared to theory (solid lines) developed using 
partially coherent optics - Equations [14] and [25] For both sets 
of data, an orange LED (A ~ 600 nm) was placed directly 
in front of a millimeter- wide slit, which illuminated a second 
slit of the same width 1 meter away. Intensity profiles were 
obtained using a CCD placed 10 cm behind the second slit. 



model (GSM) can be used to accurately predict the char- 
acteristic width, the transverse coherence length (coher- 
ence width) , and the divergence angle of a slit-collimated 
beam. In [1], we show how a GSM beam can be used 
to efficiently analyze the role of partial coherence in a 
variety of grating interferometers, and we suggest this 
model can now be applied to several matter wave inter- 
ferometry experiments. To this end, the purpose of this 
manuscript is to show that GSM beams can accurately 
model slit-collimated beams commonly used for matter 
wave experiments. 



II. THEORY - PARTIALLY COHERENT BEAMS 
FROM TWO CONSECUTIVE APERTURES 



Figure 1: In this paper, we assume that an aperture with char- 
acteristic width so (either a slit or a Gaussian filter) placed 
at z = zo is incoherently illuminated. A second similar aper- 
ture with width si is placed a distance zoi away. Generally, 
beams resulting from this arrangement are partially coherent; 
the characteristic size of the mutual coherence function (the 
coherence width) is less than or equal to the width of the 
intensity distribution. 



Here we introduce the general formalism we will use. 
We first introduce the conventional form of a function 
that describes two-point correlations in a field. A semi- 
coherent scalar field can be described as a statistical dis- 
tribution of coherent beams. One way to model this is 
to use the mutual intensity [2[, defined as 



J(pa,Pb) = (ljj(paW(pb)) 



(1) 



where the angle brackets denote an average over the en- 
semble. This expression is a way to quantify the correla- 
tion between two points, p a and p&, in a particular plane 
transverse to the direction of propagation. This can be 
seen in an alternative form for the mutual intensity: 



J(pa,Pb) = y/l(fta)\ I(Pp)v(Pa:Pb) 



(2) 



where fi(p a ,Pb) is the spectral degree of coherence, also 
defined in [2J . The intensity of the light field at a point 
p can be found by setting the two points of the mutual 
intensity equal to eachother: 



i(p) = j(p,p) 



(3) 



For use with polychromatic fields, such as those com- 
monly encountered in matter wave beams, the mutual 
intensity may also be used to calculate the spectral den- 
sity, as noted by Wolf and Devaney [3] . 

The mutual intensity provides a way to use techniques 
developed for coherent optics to examine the evolution 
of semi-coherent optical fields. The difference is that the 
semi-coherent calculation requires two operations (one 
for each of the two correlated points) whereas the conven- 
tional coherent calculation requires only one. For exam- 
ple, when calculating how a thin optical element modifies 
a mutual intensity distribution, the complex transmission 
function t{p) describing the element appears twice: 



J(p a ,p h ]Z Jr ) = t(p a )t*(p b )J(p a ,p b ;z) 



(4) 



where z + denotes a position just after the element. Like- 
wise, one can use the paraxial (Fresnel) approximation 
to the Huygens- Fresnel principle to propagate a partially 
coherent wave from a source plane zq over a distance zqi 
to plane z\ j4[: 



A. 'Hard-edge slits' 

Let us now use the above framework to calculate the 
partially coherent beam emerging from the two aperture 
system illustrated in Figure 1, assuming slits with hard 
edges. We will only calculate the mutual intenstiy in 
one transverse direction x. The complex transmission 
function for the first slit with hard edges is 



to(x) — rect I — 



(6) 



where sq is the width of the slit and, as in [6], the rect 
funtion is defined as 



rect 



, r , . 0, \x\>s/2 

(?) = h M = «/2 

W ' 1, \x\ <s/2 



(7) 



When this aperture is placed in front of a completely 
incoherent source with radiant emittance io, at zq, the 
mutual intensity directly behind the slit, at z$ , is given 
by 



J(x a , x b ; z£) = 7 rect I — j rect ( — j S(x h - x a ). (8) 

The incoherent nature of the source is reflected in the 
presence of the delta function; i.e. there is no correlation 
between any two points. 

The 1-D form of Equation [5] is used to calculate the 
field a distance zqi away. Making use of the delta function 
in Equation[8l and recognizing that rect(x) 2 = rect(x) we 
find: 



J(Pa,Pb',zi) = 
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A 2 ; 



X J(Pa>Pb', Z oWa d Pb' 



X7^[(pa-Pa) 2 -(Pb-Pb) 2 



(5) 



where p' a and p' h are points in the source plane. The 
Fresnel propagator we have chosen to use - the quadratic 
phase factor - approximates hemispherical Huygens 
wavelets with paraboloids, and it is valid for calculat- 
ing fields not just far from the source plane but also near 
to it. Equation[5] is not the van Cittert-Zernike theorem, 
which is restricted for use with fairly incoherent sources. 
It is rather a paraxial approximation to Zernike's more 
general propagation law for fields with arbitrary coher- 
ence properties [5[. 

Just as in coherent optics, the propagation of partially 
coherent fields through complex optical systems may be 
computed using succesive applications of Equations [Hand 
[51 For the remainder of this manuscript, we will only 
evaluate the beam's mutual intensity in a single trans- 
verse direction. 
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(9) 



The remaining integral over x' a becomes a straightfor- 
ward Fourier transform of a rect function and we arrive 
at: 



J(x a ,x b ;zi) = /isinc 
where I\ = IoSo/Xzoi and 

sinc(ax) 
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(11) 



We now let this mutual intensity pass through a second 
aperture, another hard-edged slit of width si: 
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(12) 



The mutual intensity just after this aperture (denoted by 
z^) becomes: 
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(13) 
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Now we wish to determine the properties of the beam 
an arbitrary distance away from the second collimating 
slit. We may again apply the paraxial approximation to 
Zernike's propagation law (Equation [5]): 
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(14) 



Even after rearranging terms to reveal a two dimensional 
Fourier transform, it is difficult to spot a full analyti- 
cal solution to Equation [TH This is unfortunate since in 
many optical and matter wave experiments great insight 
could be gained by analytically computing how such a 
beam evolves through a given system. It is for this rea- 
son that we suggest using a GSM beam to model these 
experiments. The Fourier transforms present in Equa- 
tion [TH are conducive to numerical evaluation by a com- 
puter, however, and later in the manuscript we will use 
this technique to compare this slit-collimated beam to a 
GSM beam. 



J(xa,x h ]z£) =he 7 o a b S(x b -x a ). (16) 

As before, the mutual coherence function after this first 
aperture is a delta function, indicating there is no corre- 
lation between any two points. 

Like in the previous section, we now insert Equation 
[T6l into the 1-D form of Equation [5] to calculate the field 
emitted by this Gaussian source a distance zqi away: 



J(x ai X h ]Zi) 
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xe s o " ' 5(x' b — x' a )dx' a dx' b (17) 



After making use of the delta function by integrating 
over x' h , the remaining integral becomes a straightforward 
Fourier transform of a Gaussian function and we arrive 
at: 



J(x a ,x b ;zi) = 
where h = IqSq/Xzqi, 



^riXh—Xa) i7r ( \ 



n = zoi 

is the radius of wavefront curvature, and 



cri 



Xzoi/sq 



(18) 



(19) 



(20) 



is the correlation width. We will discuss the physical 
meaning of the correlation width shortly. 

The field described by Equation [18] is incident upon a 
second Gaussian- weighted aperture with width si located 
at z\\ 



B. Gaussian apertures - GSM beams 

Again, we now calculate the partially coherent beam 
emerging from the system illustrated in Figure [TJ but this 
time assuming the two Gaussian apertures. In order to 
compare the dimensions of the resulting GSM beam to 
those of the slit-collimated beam, we must first choose 
the size of the Gaussian apertures carefully. The trans- 
mission function of the first aperture is 



t (x) 



(15) 



where we choose a convention for the width so, which is 
just slightly larger than the FWHM, such that the total 
area of the Gaussian aperture is equal the the total area 
of a hard-edged slit with the same width [6] . 

The mutual intensity directly behind this incoherently 
illuminated Gaussian aperture is given by 



h(x) 



(21) 



The mutual intensity just after this second aperture (de- 
noted by Zi) can be found using Equation [4j 



J{x a ,x h ',z J [) = ti(xa)t*(x b )J(x a ,x b ]Zi) (22) 

T ^r(xl+xl) f^I^b-X a ) 2 ™-(xl-xV) 

We now define the transverse coherence length or coher- 
ence width l\ 



(23) 



so that we may write the partially coherent beam exiting 
from the second aperture in a more recognizable form: 



J \X ai X\y\ Z-^ ) — L\& 
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Equation [2H is the usual representation of a Gaussian 
Schell-model (GSM) beam [J 0. The GSM beam can 
be thought of as an ensemble of coherent beams each 
with an average momentum and radius of curvature that 
depends on its displacement from the optical axis. 

The correlation width o\ is commonly touted as the 
characteristic maximum distance between any two points 
in plane z\ for which there is significant correlation in 
the field [8]. However, this statement loses its meaning 
when the correlation width exceeds the size of the classi- 
cal intensity distribution (a feat easily accomplished by 
choosing a sufficiently small width for the first slit or sec- 
ond slit). A hypothetical Young's experiment, which is 
the most direct method for measuring two point correla- 
tions, with slit separations much greater than the extent 
of the beam would not yield visible interference fringes. 

The coherence width i\ is more directly related to ex- 
periment, as it can never be larger than the transverse 
extent of the beam. The coherence width is the maxi- 
mum separation between slits in a Young's experiment 
for which one can expect to see large intensity modula- 
tions in the fringe pattern. Thus, we elect to use Equa- 
tion [22] as the form for the GSM beam. 

Upon another application of Equation O one may de- 
rive an analytical solution for the GSM beam at any dis- 
tance from the second aperture: 
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(26) 
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(28) 



III. 



MODELLING A SLIT-COLLIMATED BEAM 
WITH A GSM BEAM 



As stated previously, while Equation [TH provides a way 
to compute the mutual intensity of the partially coherent 
beam emerging from collimating slits exactly (within the 
limits of the Fresnel approximation) , it is difficult to cal- 
culate how this beam interacts with any optical elements 
such as diffraction gratings. In matter wave interferom- 
etry, where one is often analyzing the occurence of in- 
terference fringes, knowledge about the detailed shape of 
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Figure 3: Comparison of slit-collimated and GSM beam in- 
tensities and mutual coherence profiles. Here we model our 
own sodium atom beam, using A = 17 pm for 1000 m/s atoms, 
and 10 /im slits spaced 1 meter apart. The intensity profile 
of the slit-collimated beam (solid red) was obtained numeri- 
cally using Equation 1141 and setting x a — Xb — x. The mu- 
tual coherence profile at the center of the beam (dotted red) 
was obtained similarly, setting x a — —Xb and averaging over 
several wavelengths. The intensity profile of the GSM beam 
(solid blue) is just a Gaussian curve using Equation [26] for the 
width, and likewise for the mutual coherence profile (dotted 
blue) using Equation [27l 



the intensity distribution or mutual coherence function is 
often unnecessary. What is particularly useful, however, 
is knowledge of the characteristic size of these distribu- 
tions - the intensity and spatial coherence widths, and 
radius of wavefront curvature. The GSM beam provides 
a much faster way to calculate these properties, and a 
simpler way to analytically propagate a partially coher- 
ent beam through a given optical system. 

Comparing this field directly behind two consecutive 
slits to the field behind two consecutive Gaussian aper- 
tures (Figure |3(a)[ ), we of course see that the inten- 
sity distribution is entirely determined by the uniformly- 
illuminated second aperture. Note that the intensity of 
the GSM beam directly behind the second aperture is re- 
duced to about 20% of its maximum value at x = ±si/2, 
whereas the intensity of the slit-collimated beam is re- 
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Figure 5: Characteristic widths of the same beam profiles 
shown in Figure |4] The slashed red line displays the width of 
a hard-edged slit-collimated beam, the solid blue shows the 
width of a corresponding GSM beam, and the red dotted line 
and blue dashed line indicate the coherence width of the slit- 
collimated beam and the GSM beam, respectively (evaluated 
at the center of the beam) . 



Figure 4: (a) Intensity profile and (b) mutual coherence profile 
of a beam emerging from a hard-edged collimating slit, as a 
function of the propagation distance. Comparable profiles 
(c) and (d) are shown for a GSM beam. These figures were 
generated with the same parameters used in Figure [3] 



duced by half. When comparing profile widths for differ- 
ent types of beams we use the characteristic width of the 
profile when it is 20% of its maximum. 

The ratio of transverse coherence length to beam width 
in Equations [26] and [27] remains unchanged as the beam 
propagates through free space. That is, 



f3 = £(z)/s(z) = constant 



(29) 



where (3 is the 'degree of global coherence' @, 0|. If /? is 
equal to unity, Equations [26] to [28] describe the beam di- 
ameter and radius of wavefront curvature for a standard 
(fully- coherent) Gaussian beam. According to Figure [5] 
the global degree of coherence of our sodium atom beam 
used for inter ferometry experiments is about /3 — 0.1. 



IV. CONCLUSION 

In summary, we have shown that a Gaussian Schell- 
model (GSM) beam can efficiently model various proper- 
ties of a slit-collimated beam. While the specific shapes 
of the intensity and coherence profiles can differ dramat- 
ically between the two types of beams, the character- 
stic size of these distributions are identical and evolve 
through space via propagation in a similar way. The pri- 
mary advantage of using the GSM beam to model the 
complicated dynamics of a beam generated from hard- 
edged apertures is that it is much easier to analyze its 
propagation through various optical elements, such as 
those used in matter wave inter ferometry experiments. 
As an example of this, in an upcoming paper [1] we use a 
GSM beam to efficiently model a wide variety of grating 
interferometers commonly used in experiments, and the 
role that partial coherence and beam divergence has on 
their behavior. 
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